Solitosynthesis induced phase transitions 



Lauren Pearce 

Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1547, USA 

We consider a phase transition induced by the growth of Q-balls in a false vacuum. Such a 
transition could occur in the early universe in the case of broken supersymmetry with a metastable 
false vacuum. Small Q-balls with a negative potential energy can grow in a false vacuum by accretion 
of global charge until they reach critical size, expand, and cause a phase transition. We consider the 
growth of Q-balls from small to large, using the Bethe-Salpeter equation to describe small charge 
solitons and connecting to the growth of larger solitons for which the semiclassical approximation 
is reliable. We thus test the scenario in a simplified example inspired by supersymmetric extensions 
of the standard model. 



INTRODUCTION 

Q-balls jTj are non-topological solitons [5J [5] that are 
stable because they carry a conserved global charge. 
They arise in a number of models, and, in particular, in 
supersymmetric extensions of the Standard Model, where 
they carry baryon and/or lepton number [1]. Stable su- 
persymmetric Q-balls can form in the early universe from 
the fragmentation of Affleck-Dine condensate [S] or in 
other processes [7], and they can play the role of cosmo- 
logical dark matter [5l|8]. Furthermore, it has been sug- 
gested that Q-balls can facilitate phase transitions even 
when the tunneling rate is too small for the phase tran- 
sition to occur otherwise; the Q-balls accumulate charge 
until they reach a critical charge, at which point they ex- 
pand, causing a phase transition [9]. Such a phase tran- 
sition could have interesting cosmological implications. 

While the possibility of such a phase transition has 
been explored in the literature [TOl E] , a complete model 
of it has not yet been demonstrated. This is due to diffi- 
culties with the quantum nature of small charge Q-balls 
and also with the properties of Q-balls in the false vac- 
uum. This paper will demonstrate, from beginning to 
end, a model in which a phase transition is induced by 
solitosynthesis of Q-balls. 

This paper is organized as follows: first we specify the 
potential that gives rise to our Q-balls and show that it 
has the requisite properties; then we consider the prop- 
erties of the non-topological solitons in the false vac- 
uum. There are primarily four regimes to consider. For 
large charges, the thin wall semiclassical approximation is 
valid, while for smaller charges, the thick wall semiclas- 
sical approximation is valid. For intermediate charges, 
we interpolate between these two regimes. For extremely 
small charges, quantum effects are important and the 
semiclassical approximation is invalid; instead, we apply 
the Bethe-Salpeter equation. 

After we have described the radii and energies of the Q- 
balls, we proceed to consider the properties of the phase 
transition; in particular, the critical charge and the crit- 
ical radius. Then we consider solitosynthesis, the pro- 
cess by which Q-balls grow by accreting of charge. We 



find the temperature at which such growth begins, and 
then we calculate the rate of growth in each regime. We 
demonstrate that the growth is not hindered by charge 
depletion and freeze out, which could end solitosynthesis 
before critically sized Q-balls form. Finally, we discuss 
explicit numerical examples to show that such a phase 
transition is a theoretical possibility. 

In all this, we use a simplified toy model inspired by 
the Minimal Supersymmetric Model (MSSM). In the last 
section, we discuss the application of this analysis to the 
MSSM, and in particular we consider phase transitions 
of cosmological interest. 

THE POTENTIAL 

For this paper, we will use an MSSM-inspired poten- 
tial 

U = + "fh^ Aohf + ^fh^ + + 

(1) 

in which g is a squark field and h is the lightest Higgs bo- 
son. For simplicity, we use real fields. In general, renor- 
malization effects, including the effects of Q-balls if they 
exist, can be significant [H]; we take the couplings to be 
the renormalized couplings. We assume that any other 
particles which carry baryon number are heavier than 
the squark, to ensure the squark's stability. The origin 
is always a local minimum of this potential; however, 
for particular values of the coupling constants a different 
global minimum exists. For example, if = 200 GeV, 
ruh = 10 GeV, Aq = 240 GeV, Ai = A2 = .1, and 
A3 = 19, the origin is a local minimum and there are 
global minima at (h) — 169 GeV and (q) = ±631 GeV. 
In Fig. [1] we show a contour plot of this potential. 

Phase transitions involving multiple fields are difficult 
to solve exactly; a reasonable approximation is that they 
occur along the line connecting the two minima. This is 
a valid approximation when the potential does not have 
an unusually shaped barrier between the minima; we see 
no sign of an unusual barrier in the contour plot. The 
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FIG. 1. A contour plot of the potential U{q,h). There is a 
local minimum at the origin and a global minimum at (h) = 
169 GeV and (g) = 631 GeV. 

Potential Along Line Connecting Minima 




FIG. 2. This is the potential along the line connecting the 
false vacuum to the true vacuum. We see a large barrier, 
which suppresses phase transitions driven by thermal fluctu- 
ations. 



potential along this line may be found by substituting 
q = 0sin(6l) and h = 0cos(6') with 9 = 1.309: 

C/(0) = 18700 GeV^ c/)'^ - 58.0 GeV c/}^ + .04470"* (2) 

which is shown in Fig. [2] We see that there is a large 
barrier between the minima; this suppresses phase tran- 
sitions induced by thermal fluctuations. 

Next we must demonstrate the existence of Q-balls in 
the false vacuum (located at the origin). The conserved 
charge carried by our Q-balls will be baryon number; this 
is conserved in the false vacuum because squarks do not 
have a vacuum expectation value. The condition for the 



existence of Q-balls involving multiple fields is that 211/ 
J2k Qk4''i is minimized at a nonzero value of the fields, 
where Qt is the charge of of the field 0fe 1110]. Because 
the baryon number carried by the squark field is 1/3, we 
consider the minimum of 611 jq^ . Due to the q^h term, the 
origin will not be even a local minimum of U /if'] in fact, 
the global minimum is located at go = 624 GeV and — 
168 GeV, and in which the potential is -6.18 • lO'^ GeV*. 

Therefore, Q-balls carrying baryon number exist in 
this vacuum, and furthermore, the potential is negative 
within their interior. This is necessary for a Q-ball in- 
duced phase transition to occur; such a phase transition 
converts the fields everywhere to their values inside of 
the Q-ball, qo and h^. Since this is not quite the true 
minimum of the potential, the system will then slide clas- 
sically into the true minimum. 

Again we consider the potential along the line connect- 
ing the initial vacuum and the final state, which, with 
q ~ 0sin(6'), h = (pcos{0), and 9 = 1.308, is: 

J7((/)) = 18700 GeV^ 4>^ ~ 58.2 GeV (^^ + .04500^ (3) 

We can also describe the Q-balls in terms of the field. 
The field inside such a Q-ball is the value of 4> that min- 
imizes U{4))/(t)^, which is 00 = 646 GeV, and as ex- 
pected, this is Y^^o + ^0- Also, as expected, C/((/)o) — 
—6.18 • 10^ GeV*. Thus, at least for large charge, we can 
consider our Q-balls as being bound states of squarks ex- 
changing Higgs bosons, or as coherent oscillations of 
quantum. For small charge states, we should remember 
that the charge is really carried by squarks that exchange 
Higgs bosons. 

Finally, because squarks carry charge 1/3, the charge 
present in a field configuration is: 

^ uj f ^3 .2 sin2(1.308) r 3 ^2 /.N 

where w describes the oscillatory time dependence of the 
fields. For later convenience, we measure the charge in 
units such that Q = sin^(1.308)Q'/3; then: 

Q' = Lo j d^xcj)^. (5) 

Physically, Q' is the charge in a single quanta. In 
terms of units of charge Q', a single squark carries charge 
3.107; this will be even closer to 3 the nearer the global 
minimum is to the q axis. Speaking loosely, since 9 ^ tt/ 
2, the field (p is "almost" the squark field, and the charge 
Q' is approximately the number of squarks present. 
Generically, then, we consider a potential of the form 

c/(0) = ^0'-^0'+^0^ (6) 

in which the field is made of squarks and higgs bosons, 
and furthermore, the charge of the field is approxi- 
mately the charge of the squark. The general form of the 
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relation between the charges is: 

3Q 



sin' (9) 



(7) 



Corrections to the Potential 

At the scale of color confinement, we expect the 
squarks to arrange into color singlets of the form 
^abc^apQaQb^c^ where Greek letters denote SU{2) in- 
dices and Latin letters denote color indices [T3]. To 
avoid this complication, we will ensure that the relevant 
temperatures are above the scale of the QCD phase tran- 
sition; furthermore, the scalar binding interaction medi- 
ated by Higgs bosons is much stronger than the strong 
interaction. 

In general, finite temperature corrections to the po- 
tential, \^T^(fP' and XTcj)^, may be important; however, 
when the relevant temperatures are significantly smaller 
than mo and A, these finite temperature effects can be 
neglected. This is the case in the examples considered. 



PROPERTIES OF Q-BALLS IN THE FALSE 
VACUUM 



write: 
E{Q') 



Q 



12 



Ui(f>) d^x 



-T + V, 



(11) 
(12) 



where T ^ j {XJcjif (fx/2 and F = / U{(j)) (fix. 

In thin wall approximation, w for r < R — 5/ 
2 and w for r > R + 5/2, where 5 is the width of 
the surface of the Q-ball. The thin wall approximation 
is valid if (5 <C i?. The "volume" term then has two 
contributions, one from the interior of the Q-ball and 
one from the surface. In the interior of the Q-ball, the 
potential is U{(j>o), while in the surface of the Q-ball, it 
is /3mo0g, where /3 is a positive constant. Therefore the 
"volume" term is: 

[ U{4>) (fx = ^7r[/(0o)^^ + 47r(5i?2 • /3mg0o- (13) 



In the surface, the field changes by A(/) — (j>o in the 
distance Ar = S, thus d(j)/dr <j)o/5. Introducing a 
constant a to account for the uncertainty, the surface 
term is: 



The exact and general equation for the energy of a 
Q-ball of arbitrary charge has three terms: 



E(,Q') 



d'x{\\^\' + \\Vcl^\^ + U{c^)]. (8) 



The field oscillates in time as e^'^^(j){x) where the fre- 
quency is related to the charge according to: 



[d^x(i>*'dt 

2t 



oj I (f)'d^x 



(9) 



After some manipulation, one can write |16j : 



E 



(10) 



where the first term is the three-dimensional Euclidean 
action of the "bounce" solution tunneling between the 
two minima of the effective potential Uui{4') — U{4)) — 
c^20V2. 



0^ 



(14) 



The energy must be a minimum with respect to both b 
and i?; minimizing with respect to 5 gives 8 — y/aj]3m^ 
and 



E 



— ^+8^mov^-i?20^-|--^t/(0o)i?^ (15) 



By manipulating equation (10), one can relate ^/aP to 



the one-dimensional Euclidean action for the true poten- 
tial: 



Si - 



y^W{^ d(j) = 2mo^/^ (j)l, (16) 



which is related to the three-dimesional action through 
53 =47ri?2c/(0o)/3 + 47ri?2S'i [ig. 

Minimizing the energy with respect to R results in a 
constraint between the charge and the radius: 
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+ 167rmo + 47rC/((/)o)i?^ (17) 



Thin Wall Regime 

In the thin wall regime, the energy may be calculated 
in a different manner, due to [14) . Beginning again with 
equation ([s]), we use the oscillatory time dependence to 



If U (0o) > 0, the last term would dominate over the sec- 
ond term, and if we neglect the second term, we could 
solve for R in terms of the charge. This would give 
the familiar R cx Q'^/^ behavior [T]. However, we have 
L/(0o) < 0, and so one cannot neglect the second term. 
This sixth order equation has no closed form solution. 
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Thick Wall Regime 



Bethe-Salpeter Regime 



For the thick wall approximation, we consider equation 



(10 1. As the charge becomes small, the frequency ui be- 



comes large; then the true asymmetric minimum in Ui^ (0) 
is significantly lower than the symmetric minimum; this 
is true independent of the sign of C/(0o); in fact, since we 
have U{4>q) < 0, the second minimum (the true vacuum) 
is lower than the symmetric minimum even at w = 0. 
Therefore, the Euclidean action is the same in both 
cases, and the relations between the energy, radius, and 
charge are unchanged. These are |16j : 



E = Q'mo 1 



6 8 



(18) 
(19) 



where e = Q'A^/SS'^mg and S^p 
numerically. This is valid when: 



4.85 was determined 



35'^TOo 



Q' < 



2^2 ■ 



(20) 



The thick wall approximation, like the thin wall 
approximation, neglects quantum corrections; thus it 
breaks down when these are large, which occurs around 
Q' < 7 [T7]. 



Intermediate Regime 

There is a regime between where the thin wall approx- 
imation is valid and where the thick wall approximation 
is valid; unfortunately no general solutions are known in 
this regime. Therefore, we use a linear interpolation be- 
tween the two regimes. We will need only the radius in 
terms of the charge Q' , for which we use: 



R 



Q'~7 
Q' 



Rt 



7 



hin 



:R 



■thick ■ 



(21) 



We have already written Rthick in terms of Q'; how- 
ever, as found above there is no closed-form equation for 
Rthin in terms of Q' . Therefore, we use a numerical ap- 
proximation for the thin wall radius at small charge of 
the form i? « a + bQ' ^Z^. While this can only be jus- 
tified by its numerical accuracy, it can be motivated by 



neglecting the third term in the constraint equation (17 1. 



Then in the intermediate regime the linear interpolation 
becomes: 



R ^ 



Q'~7 
Q' 



bQ 



I 2/5 



215', 



4, 



Q' ■ 



(22) 



The stability of states with very low charge is vital to 
building critically sized Q-balls, but for these states, one 
cannot use the approximations already discussed because 
quantum effects are important. At small charge, it is 
furthermore incorrect to think of the Q-balls as made of 
quanta; instead we should remember that they are bound 
states made of squarks and Higgs bosons. We will do this 
until the bound states are large enough that quantum 
effects are negligible, which occurs at Q' = 7 [T7] . 

For these small Q-balls, we consider them as bound 
states of squarks exchanging light Higgs bosons; this is 
approximately described by the Bethe-Salpeter equation 
in the ladder approximation. This neglects diagrams 
where the rungs of the ladder are crossed; therefore, we 
use an effective coupling A which we tune to ensure that 
the energy at large charge matches the result from the 
thick wall approximation. The Bethe-Salpeter equation 
is discussed in more detail in the Appendix. 

In this section, we work with the number n — SQ of 
squarks present in the state; at the end, we will relate this 
to the charge Q' that we have been using by Q' = 3Q/ 
sin^(0), as in ([t]). The lowest state is a single squark, 
and then the first step is the relatively simple case of two 
equal squarks forming one n — 2-ball; in this equal mass 
and equal coupling case the Bethe-Salpeter equation (af- 
ter a Wick rotation) is: 



M 
Y 



M 

i2 

16^ 



ip{p) 



{p - A:)2 + ' 



(23) 



where tpip) is the wavefunction. Approximating ~ 0, 
the bound state energies are [18j : 



M„ = 2m 1 - 



8n2 



= 2m„ 



1 - 



A* 



20487r2m4n2 



(24) 

if a = /WTrrrig < 1. We may find the binding energy 
for the ground state by taking n — I. 

For the remaining states, the masses and couplings at 
the top and bottom of the ladder are unequal; after a 
Wick rotation, the Bethe-Salpeter equation in the ladder 
approximation is: 



[{m + Af + {p^ 



irjiPf] [(m- A) 
_ A^ 



[p + iri2Pf] 



(p ~ ky 



2p{p) 

(25) 



where the masses of the particles, rritop and rribottorm 
are m ± A. The coupling constant is A' = gtopgbottom/ 
16n^ ^ (n — 1)A2/I67r2, where gtop = (n — 1)A is the cou- 
pling at the top of the ladder and gbottom = ^ is the cou- 
pling at the bottom of the ladder. The total charge Q of 
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n (No. Squarks) 


Mass of Q-ball 


Binding Energy 


1 


m. 





2 


2m, - .0000989i''/m^ 


.0000989i''/m^ 


3 


Sm, - .0002309i''/m|| 


.0001320i''/m^ 


4 


4mg - .0003793i''/m^ 


.0001484i''/m^ 


5 


5mg - .0005376i''/m^ 


.0001583i''/m^ 


6 


6m, - .0007025i''/m^ 


.0001659i^/m^ 


7 


7m5 - .00087211^/m^ 


.0001696A''/m^ 



TABLE I. Energies of small Q-balls from the Bethe-Salpeter 
equation. We note that the dependence of the mass on A 
matches the dependence of the mass in the thick wall regime 
on A, and similarly the dependence on m, matches the de- 
pendence in the thick wall regime on mo. 



Comparing the second terms in the mass equations 
gives A = 4.51A. Using the value of A in the potential in 
([3]) gives A = 263 GeV, while the value that we put into 
our original potential in ([T]) is 240 GeV. This is a signif- 
icant difference, which we attribute to a combination of 
the inherent inaccuracy of the ladder approximation in 
the Bethe-Salpeter equation; it does not improve if we 
iterate further. 

The equation for the masses given above is valid as 
long as A' <C m^. For n = 7, this is approximately 7A'^ / 
167r^ • ArUq — .019, and so the approximation remains 
valid. 

CRITICAL VALUES FOR PHASE TRANSITION 



the resulting Q-ball is n/3. The energy-momentum four- 
vector of the bound state, P, is given by (0, M) where M 
is the bound state mass. 771 and 772 come from transform- 
ing to a "center of momentum" reference frame; these 
are: 



mto 



^top H~ '^bottom 



m 



'^bottom 



The binding energies are [19]: 



4A^ + 4m^ 1 



A2 



Tfltop + rnbottom 



4m4 (l-A2/m2)2 
(26) 



as is explained in more detail in the Appendix. We use 
this equation iteratively to find the masses and binding 
energies of the small charge Q-balls; the results are shown 
in Table U 

We have calculated this until n = 7, or Q = 7/3. Using 
6 — 1.308, this corresponds to Q' — 7.51. (Recall that 
we are measuring charge in units of the charge of 4> and 
each squark has slightly more than 3 unit charges.) Since 
this is greater than 7, the thick wall approximation is 
applicable [T7j. In the thick wall regime, a Q-ball with 
this charge has energy 

M = 7.51mo - .333A*/ml. 

The difference in the first terms comes from the fact that 
states described by the Bethe-Salpeter equation do not 
have exactly the same proportion of squarks and Higgs 
bosons as the Q-balls described in the thick wall regimes; 
however, the difference in these terms is 7.2 percent. This 
can be further improved by moving the global minimum 
closer to the q axis. This is because: 

Q'mo = n y'TO?sin2(e') +m2 cos2(6l)/ sin2 (9) , (27) 

as can be seen by comparing the potentials ([T]) and (|3]), 
and using the relation between n and Q' . 



In the thin wall regime, the interior of the Q-ball is 
in the true vacuum, which has negative energy density. 
If charge continues to increase, the Q-ball expands, con- 
verting more of the space into the true vacuum. At a 
particular value of the charge and radius, it expands un- 
controllably, thereby converting all space into the true 
vacuum _9J. This Q-ball induced phase transition can 
occur even when such a phase transition cannot be in- 
duced by thermal fluctuations. 

As will be demonstrated in our numerical example, the 
critical charge is of order 10^, which is within the thin 
wall regime. At the critical point, not only is dE/dR = 
0, but also d^E/dR^ = 0, which gives the additional 
constraint 



= 



2^4,1 



+ lQTTmQ^Jal3(t)iRl ~ SttUoR^, 



(28) 



where U{(t>Q) = —Uq with Uq > 0. We solve the two 
constraint equations (17) and (28) for the critical charge 
and critical radius: 



Rr = 



10mo-\/a/ 



3C/n 



t = 27rfi 



8U0 /IQmoVaP^ 
45 I 3C/o 



(29) 
(30) 



SOLITOSYNTHESIS 



In thermal equilibrium, the number density of Q-balls 
of a particular charge is given by a Saha equation: 



94>9Q' 



\moTj 



(31) 



where Bqi is the binding energy of a soliton of charge Q' 
and gqi is the internal partition function of the soliton. 
50 is the number of degrees of freedom associated with 
the (j) field. Noting that we chose to work with real fields, 
this is 3 from the color charge carried by (j). n^, also called 



6 



the charge density, is the number of free squarks; since 
9 « 7r/2, conceptually we can think of this as the number 
Q' = 1-balls present. This is given by: 



710 = ririj 



(32) 



'>2 



where the baryon asymmetry is 77 and in any radiation- 
dominated era the photon density is 2C(3)r'^/7r^ |20j . 

The typical approach would be to solve these coupled 
equations numerically. However, the critical charge is 
generically of order 10'^ to 10^, which leads to at least 
10^ coupled equations. It is infeasible to solve these si- 
multaneously. Therefore, we take a different approach 
following [5T] and consider the evolution of the single Q- 
ball. One significant advantage is that we will see that 
the Q-ball grows fast enough that we can ignore charge 
depletion and set 



2.404r3 



(33) 



It is possible for the binding energy to be sufficiently 
large that solving this equation for temperature results 



in an imaginary value; returning to equation (36), this 



happens when Bq is so large that the right hand side 
is always less than 1, which means that ngi > tiq'-i 
always. Physically, whenever a {Q' — l)-ball forms, it 
will always grow into a larger ball; we may say that the 
solitosynthesis temperature for these charges is infinite. 

Smaller charges have smaller binding energies, and so 
we will find real solitosynthesis temperatures for small 
Q-balls. Therefore, we will need to wait for these smaller 
Q-balls to form, and then wait for these to grow into the 
larger ones which can always grow. 

As this suggests, Tg is greater for larger charges, which 
we will verify numerically. Therefore Q-ball growth is 
a winner-takc-all-situation, and the solitosynthesis tem- 
perature cannot cut off a growing Q-ball. 



Rate of Diffusion 



A single Q-ball grows or shrinks according to: 



dt 



'^absi^Q ) ^evapiQ )i 



' evap 



(34) 
is the 



where Vabs is the absorption rate and 
evaporation rate. By detailed balance, nqirabsiQ') — 
UQi+irevapiQ' + 1); also, the rate of absorption is Vabs — 
n^v^UabsiQ')- For large charges, Uabs ~ t^R^- We will 
see numerically that the radius does not change rapidly 
as a function of charge; then (Jabs{Q') ~ <^abs{Q' — 1)- 
These approximations give: 

dt 



n^V^UabsiQ') 1 - 



(35) 



Thus the determining factor is uq' -xjuQi: if it is less 
than one, absorption dominates, but if it is greater 
than one, evaporation dominates. From the Saha equa- 



tions (311, this important ratio is: 



507 



(1-kT 



-3/2 



(36) 



nqi ' ' 2.40477 \ tt^q 

At large temperatures, the exponential is negligible and 
this scales as T^^^^. However, this ratio is less than one 
only if T > r]~'^^^mQ, which is typically quite large, above 
the temperatures at which the supersymmetry which in- 
spired our potential is typically broken. Therefore, we ex- 
pect evaporation to dominate at the temperatures when 
the Q-balls are formed. 

As the temperature decreases, the exponential term 
is no longer negligible. Because —Bqi < 0, this term 
decreases ngi^i/nQi . Therefore, at some temperature Tg 
absorption will dominate. The ratio nqi-i/nQi is equal 
to one at: 



Hg^) - ln(77) -t- (3/2) ln(mo/r,) - 1.34' 



(37) 



A Q-ball grows by absorbing the nearby charge. If the 
charge is not be replenished sufficiently quickly through 
diffusion, there may be a local depletion of charge near 
the Q-ball which limits its growth. If this occurs, the 
rate of growth will be given by '"di// , the rate that free 
squarks diffuse into the surface of the Q-ball, instead of 

'^abs • 

Reference [22] is concerned with the related process of 
the diffusion of evaporating squarks away from a Q-ball. 
The particle flux through the Q-ball surface is given by: 



where D 



dQ[ 
dt 



eg 



(38) 



4 for relativistic squarks, and 



fc w 1 was determined numerically. We need to adjust 
this equation because we are concerned with particles 
diffusing towards the Q-ball; the rate has the opposite 
sign and we multiply this by the velocity of the squarks 
because they are moving non-relativistically. Thus: 



rdiff = Vit^lGnRT 



-1 en 



(39) 



We wish to compare this to Tabs j the rate of absorption 
as approximated above; the ratio is: rdiff/vabs = 4r~^/ 
R. Perhaps surprisingly, this is small for high temper- 
atures and large for low temperatures. This occurs be- 
cause the rate of diffusion is propotional to T^/^ while the 
rate of absorption is proportional to T^/^. Even though 
diffusion is decreasing as the temperature decreases, the 
rate of absorption drops faster; therefore, diffusion will 
limit the growth of Q-balls for temperatures above 4/ 
R. For radii of order .01 inverse GeV to .1 inverse GcV, 
this temperature is of order 40 GeV to 400 GeV, which 
is significantly above the solitosynthesis temperatures. 
Therefore, diffusion will replenish the charge sufficiently 
quickly at the relevant temperatures. 
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This, combined with the winner-take-all behavior, 
demonstrates that global depletion of charge is not an 
issue, provided that most of the charge is in free squarks 
during solitosynthesis. We will verify this numerically. 



Rate of Growth in the Thin Wall Regime 

Next we consider the rate of growth of the Q-balls in 
the various regimes. For temperatures below the solitsyn- 
thesis temperature, the rate of evaporation is small, and 
we may approximate dQ' /dt — n^v^UabsiQ') from equa- 
tion (35). Since charge depletion is negligible, we may 
The Q' 



have Af = 6. Then we have: 



assume = -qn^. ine C^' — 1-balls being absorbed are 
in thermal equilibrium at T ^ toq with average velocity 
= \/lT j-KTriQ. Additionally, we use the geometric area 
7ri?^ for the cross section. In the radiation-dominated era, 
the temperature and the time are not independent; they 
are related by |20j : 



(40) 



where TV is the effective number of degrees of freedom 
of the particles in thermal equilibrium, with fermionic 
degrees of freedom weighted by 7/8. In our toy model 
with only the Higgs boson, squarks, and photons, we 



= -1.34 • 10'* GeVdT/T 



Thus our differential equation is: 



1 



dQ' 



1.34 • 1018 GeV dT 



2.4 / 2T 



(41) 



(42) 



using = 2AT^/Tr'^. 

The right-hand side involves the radius which is not 
independent of the charge; however, in the thin wall ap- 
proximation the radius cannot be written in terms of the 
charge because of the form of the 6th order equation re- 
lating them. Fortunately, one can write the charge in 
terms of the radius, and then we consider the rate of the 
growth of the radius of the Q-ball until it reaches the 
critical radius: 



Then the differential equation is: 



(43) 



5a^+6aa£ ^ _ 8.17 ■ 10^^ GeV, 

(44) 

Both sides of this equation can be integrated explicitly: 



a5 , / fl5 + 2a6i?^ + 2 Vaei?. (as + a^R,) \ 2 8.17 ■ lO^^?? GeV ^^3/2 g/aX 

I ) ~ 3 V^^*-* " ^ J ' ^^^^ 



where Ri is the radius of the smallest Q-ball at which 
the thin wall approximation is valid and T^tart is the 
temperature at which this Q-ball starts to grow. This 
can be less than Tg if these Q-balls do not form until 
a lower temperature. If we set R — Rc, this equation 
can be solved for the temperature at which the Q-ball 
becomes critically sized. 



differential equation becomes: 

§; = -7.36.10-GeV^5^Vr, (46) 



whose solution is: 



Rate of Growth in the Thick Wall Regime 



We begin with the differential equation (42) which is 
valid in the thick wall regime also. We directly relate 
the radius to the charge, R = SSti^mo/Q' A^; then the 



Q^3-Q^ 3 = 14.7-101** GeV 



3/2 
V 



^ start ^ 



(47) 

where Tgtart is the starting temperature for thick wall 
growth. This is either the Q' = 7 solitosynthesis tem- 
perature, or the temperature at which Q' ~ 7-balls form, 
whichever is smaller. 
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Rate of Growth in the Intermediate Regime 



We again begin with the differential equation (42) 



and use the linear interpolation for the radius, equation 

Q'^dQ' 



(22 1, which gives: 

Q'f 

7 



{{Q' - 7)(a + bQ' 2/5) + 21S^/Q' ^A^)' 
5.45 • IQi^ GeV?? 



,3/2 



t: 



3/2 



/mo 



(48) 



The left hand side of this equation must be integrated 
numerically. 



Bethe-Salpter Equation Regime 

We next consider the growth of very small Q-balls; in 
this regime, cross sections cannot be approximated by 



the geometrical area and so equation ( 42 1 is not valid 



However, because these are the first steps of solitosynthe- 
sis, all of the charge will be in these lowest seven states, 
which we label with n, the number of squarks present in 
the state. Therefore, one can return to the initial method 
of considering the evolution of the number densities as a 
function of temperature; we have 8 equations to solve 
numerically, instead of 10^. 

The number densities of the Q-balls are given by the 



Saha equations like (31 1, which we write in terms of frac- 



tional densities = n„n/N, where N is the total num- 
ber of squarks, 77 • 2.404X^/7:^: 

X„.^?S2f?^I)"\._,X,e=./- ,49) 



1 37r2 



mo 



with the additional equation Xi + X2 + X3 + X4 + X^ + 
Xq + Xr = 1. (The 3 in the denominator of the Saha 
equations comes from the 3 color degrees of freedom for 
a real squark field.) 

We observe that we do not need a generic Q-ball to 
grow into a critically-sized Q-ball to induce the phase 
transition, but only one per Hubble volume, 1 / where 
H = rV2.43 • lO^s GeV. Therefore, our generic ap- 
proach is to find temperature at which there are of order 
10 n = 7-balls per Hubble volume, which must be done 
numerically. Since Q' = 7.25 when n = 7, we may begin 
the thick wall analysis once sufficiently many n = 7-balls 
form. We note that we must also verify that at this tem- 
perature most of the charge remains in n = 1-balls; oth- 
erwise our analysis above is invalid because we neglected 
charge depletion. 

We see that the exponential suggests that X^ will be- 
come large at low temperature, and due to the T^^^, we 
also expect Xj to be large at higher temperatures. Thus, 
we generically expect X^ to grow at both large and small 
temperatures, with a minimum between. Typically, the 



solitosynthesis temperature for Q' = 7-balls is after the 
number of n = 7-balls per Hubble volume has dropped 
beneath 1; that is, most of the n — 7-balls have evapo- 
rated away. Then we need to wait until n = 7-balls form 
again at lower temperatures. 

However, there exist cases in which the number of 
n = 7-balls per Hubble volume is still greater than 1 when 
the universe cools to the solitosynthesis temperature of 
Q' = 7-balls; in this case, they may begin to grow im- 
mediately. In fact, then our analysis underestimates the 
temperature at which the phase transition occurs. Since 
larger Q-balls begin accreting charge earlier, it is likely 
that there are even larger Q-balls that have not evapo- 
rated away when the temperature reaches their (higher) 
solitosynthesis temperature. Again, though, we reiter- 
ate that our goal is simply to investigate the theoretial 
existence of such a phase transition. If one of these n = 7- 
balls has time to induce a phase transition, then we can 
be certain that any larger ones that had the opportu- 
nity to grow earlier would also induce a phase transition, 
and thus we still conclude that the phase transition does 
occur. 

It should be noted that especially in such a case it is 
important to verify that most of the charge is in n=:l- 
balls before we can use the results derived above; this 



means solving equations ( 49 ) numerically. 



FREEZE OUT 

Q-ball growth can be ended in one of two ways: ei- 
ther the necessary reactions freeze out as the universe 
expands, or the Q-balls deplete the nearby charge. We 
have already demonstrated that charge depletion does 
not hinder the growth of at least one critically sized Q- 
balls per Hubble volume; therefore we need only to con- 
sider freeze out. 

The reactions responsible for Q-ball growth freeze out 
when their time scale is greater than the Hubble time 
scale, th = H~^. While the universe is radiation domi- 
nated, the Hubble constant T'^/Mpi, and so the Hubble 
time scale is th = 2.43 • lO^^ GeV/T^. 

The time scale of Q-ball growth is Tabs — 1 /fats — ^/ 
n^fjv^. We consider the later reactions in the sequence; 
then the heavy Q-balls are effectively at rest and the 
Q' = 1-balls are moving non-relativistically in thermal 
equilibrium, with = [2/ ^/TT)yj2T /rriQ. We use the 
geometric cross section, u = Tri?^; then 



'^abs 



(50) 



A.SrjT^R'^ V Too ' 
Setting these timescales equal and solving for T gives 



T = 



1 



2.43 • IQis GeV • 4.877 R'^V^> 



(51) 
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Charge Q' 


Radius (GeV~^) 


Bg' (GeV) 


T, (GeV) 


Qc = 5.36 • 10^ 


.341 


9.31 • lO'' 


oo 


5000 


.0375 


6.20 • 10^ 


CXD 


1000 


.0195 


9.60 • 10* 


16600 


500 


.0148 


4.05 • 10'' 


5420 


200 


.0102 


1.16 • 10" 


1190 


7 


.118 


.482 


.0163 



TABLE II. Solitosynthesis temperatures for several charge 
values. By infinity, we mean that such a Q-ball always grows. 
The top three were calculated in the thin wall regime, while 
the last one was calculated in the thick wall approximation. 
The other two are technically in the intermediate regime; to 
approximate their solitosynthesis temperatures we used the 
thin wall regime. The important point is that because the 
temperature rise, they cannot cutoff a growing Q-ball. 

When we do our numerical analysis below, we will find 
that these are orders of magnitude smaller than the tem- 
peratures relevant to the phase transition. 



are of order 10~^ n = 7-balIs per Hubble volume. There- 
fore, we need to wait until these small Q-balls form again 
at low temperatures before thick wall growth can begin. 
Numerically, we find that there are order 10 n = 7-balls 
at T = .00889 GeV. 

Next we must address charge depletion. At the start- 
ing temperature of .00889 GeV, over 99.9999999 percent 
of the charge is in individual squarks, and so we are justi- 
fied in ignoring charge depletion. However, we also need 
to know what is the lowest temperature for which this 
assumption is valid; if our final temperature is beneath 
this, then our analysis is untrustworthy. To determine 
this bound, we consider at what temperature the major- 
ity of the charge is no longer in n = 1-balls (which are 
individual squarks), if we ignore all of the states above 
n = 7. These charge densities are shown in Figure [3] We 
see that the majority of the charge is no longer in individ- 
ual squarks around T = .0019 GeV. As long as our final 
temperature is above this, we are justified in neglecting 
charge depletion. 



NUMERICAL EXAMPLE 



Thick Wall and Intermediate Growth 



Finally, we demonstrate that the potential above is one 
in which all of these processes work out. As an reminder, 
the numbers above give mg = 193 GeV, A = 58.2 GeV, 
and A = .0450, which gives a potential where the thin 
wall approximation is valid for large charge. Above we 
also found that the minimizing field is — 646 GeV at 
which the potential is -Vq with f/o = 6.18 • 10^ GeV''. 

The other constant that must be set is 77, the baryon 
asymmetry. In the actual universe, this is about 5- 10^ 
Again, though, we emphasize that our goal is to demon- 
strate that this phase transition is a theoretical possi- 
bility, and not necessarily part of the evolution of the 
universe. Therefore, in this first numerical example, we 
will take 77 = 3 • 10~®. This value illustrates the phase 
transition well, although the final temperature will dip 
below the QCD confinement temperature; therefore we 
will also give a second, although more complicated, nu- 
merical example which avoids this. 

As regards the phase transition, the critical charge is 
5.36 • 10^ and the critical radius is .341 GeV"^ We 
present a table of the radii and solitosynthesis tempera- 
tures for various charges in Table [TTl With these radii and 
our chosen value of ry, the freeze-out temperature is of the 
order 10~^^ GeV, which is significantly smaller than any 
of the temperatures we will consider. 



Bethe-Salpeter Growth 

We begin with solitosynthesis in the Bethe-Salpeter 
regime. At the Q' = 7 solitosynthesis temperature, there 



The thick wall approximation is valid until Q' 80; 



the two constraints (20) give Q' < 80 and Q' < 228. 



The growth in this regime is virtually instantaneous; the 
temperature drops by less than one part in 10^. Thus, 
the starting temperature is still .00889 GeV for growth 
in the intermediate regime. 

The thin wall regime becomes applicable for T ^ 228; 
therefore we will use the intermediate regime for charges 
between 80 and 1000. First, however, we must find 
the constants a and b. We numerically fit the function 
a + bQ' ^/'^ to the radius for small values of Q'; the result 
is plotted in Figure m This fit gives a = —9.21 • 10 ^ and 
b = 1.24-10^'^. Substituting this into differential equation 
(48) and solving numerically gives us Tf = .00855 GeV 
for the temperature what the Q-ball reaches Q' = 1000. 
This is, of course, less than the solitosynthesis tempera- 
ture for such a Q-ball, so thin wall growth begins imme- 
diately. 



Thin Wall Growth 

First, we show the energy and radius for the thin wall 
approximation; these are shown in Figure [5j When we 
derived the solitosynthesis temperature, we assumed that 
Rq+1 ~ Rq'j we see that throughout the thin wall regime 
this is justified. Using T = .00855 GeV as the starting 



temperature in equation ( 45 ) , we find that the Q-ball 
grows to critical size at T = .00493 GeV. This temper- 
ature is greater than .0019 GeV, and furthermore, it is 
greater than the freeze-out temperature scale; therefore 
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FIG. 3. The growth of very small Q-balls. is the charge density of each type of Q-ball, given by Eq ( |49[ ). Notice that as 
temperature decreases, the number of n = 7-balls increases, while the number of n = 1-balls decreases, as we expect. 



Radius vs. Charge 
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FIG. 4. The numerical fit for the radius as a function of charge 
Q' in the thin wall approximation, for small charges. This is 
the fit to be used in the interpolation for the intermediate 
regime. 



we conclude that such phase transitions are indeed pos- 
sible. 



A SECOND NUMERICAL EXAMPLE 

As mentioned, the previous numerical example is not 
strictly speaking acceptable, because the final tempera- 
ture is less than the QCD scale, A = .217 GeV, at which 
confinement introduces additional complications. There- 



fore, we consider a second example, which avoids this 
problem. We use the same potential as above; thus mo, 
A, and A are unchanged. The critical charge and radius 
are also unchanged, as are the ranges where the thick 
and thin wall approximations are applicable. Similarly, 
the numerical fit for the radius at small charges used in 
the intermediate regime is unchanged. 

We will, however, choose an exceptionally large 77, of 
the size of 3 • 10'^. Then the solitosynthesis temperature 
for Q' = 7-balls is .258 GeV, above the QCD scale. Fur- 
thermore, at this temperature, there are still of order 
lO^'' n — 7-balls per Hubble volume; thick wall growth 
may begin immediately. This is one of the exceptional 
cases mentioned in the discussion on growth in the Bethe- 
Salpeter regime; even though these Q-balls are evaporat- 
ing away, there are sufficiently many of them for growth 
to begin immediately. For clarity, we plot the number 
density of n = 7-balls per Hubble volume in Fig. |6] to 
show it has the expected behavior. Furthermore, at this 
temperature, 71 percent of the charge is in individual 
squarks, and so we may ignore the charge depletion if 
the phase transition occurs sufficiently rapidly. 

We will show that this Q-ball can grow into a criti- 
cally charged Q-ball and induce the desired phase tran- 
sition. However, as discussed in the section on growth 
the Bethe-Salpeter regime, this may underestimate the 
temperature as which the phase transition occurs, be- 
cause larger Q-balls may have evaporated away suffi- 
ciently slowly that there are still sufficiently many of 
them to induce a phase transition when they begin ac- 
creting charge at their larger solitosynthesis temperature. 
Since there are so many n = 7-balls at their solitosynthe- 
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FIG. 5. The energy and radius as a function of charge in the 
thin wall regime, calculated numerically. 



sis temperature, this is extremely likely. However, if these 
n — 7-balls can induce a phase transition, then we can be 
certain that any larger ones that began growing earlier 
would also induce the phase transition, and so the phase 
transition certainly occurs. 

As above, the thick wall growth is virtually instanta- 
neous; there is no appreciable drop in the temperature. 
Similarly, the growth in the intermediate regime is also 
extremely fast; the temperature drops less than 1 part 
in 10^. Finally, growth in the thin wall approximation is 
equally fast; again the temperature changes by less than 
one part in 10^. Thus, such a Q-ball becomes critically 
charged within a temperature change of .000000001 GeV, 
and so the final temperature is still above the QCD scale. 

However, we can merely conclude that the phase tran- 
sition does occur; we cannot conclude that this is the 
temperature it occurs at, as we could above. As men- 
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FIG. 6. The number of n = 7-balls per Hubble volume. As 
expected, the number rises at both high and low tempera- 
tures. At r = .258 GeV, there are sufficiently many for thick 
wall growth to begin, even though they are still evaporating 
away into small Q-balls. 

tioned, it is very likely that it in fact occurs at a larger 
temperature. This is why we gave the first numerical 
example, which is a more typical case. 

POTENTIAL MSSM APPLICATIONS 

While the theoretical possibility of such a phase tran- 
sition is in itself intersting, one would also like to know 
whether such a phase transition could occur in extensions 
of the Standard Model such as supersymmetry, which 
naturally provides squarks carrying baryon number. This 
analysis suggests that, provided that the requisite vac- 
uum structure can be found, such phase transitions are 
indeed possible. 

Indeed, one can ask whether such a phase transition 
is possible within the evolution of our own universe. If 
squarks do exist, they must be significantly heavier than 
quarks. Therefore they will decay rapidly, and thus we 
cannot build critically charged Q-balls out of them in our 
current vacuum. However, such a phase transition could 
have occurred in the past, if the vacuum structure has 
these requisite properties: 

1. A global minimum in which no squarks or sleptons 
develop vacuum expectation values, so that baryon 
number and lepton number are conserved. 

2. A local minimum in which no squarks develop vac- 
uum expectation values, so that baryon number is 
conserved. 

3. In the local minimum, quarks must be heavier than 
squarks, so that they are stable against decay into 
quarks. 
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4. In the local minimum, one of the bosons that medi- 
ates an interaction between squarks must be lighter 
than the squarks; this is required for bound states 
to develop. 

5. The potential expanded in the local minimum must 
allow the creation of Q-balls through the squark 
fields. 



6. The barrier between the local minimum and the 
global minimum must be sufficiently large to sup- 
press tunneling between the minima by thermal 
fluctuations. 

Such a vacuum can indeed by found; as an example, 
consider an MSSM potential of the form: 



U 



-mlH*H + m'^^Q*Q + mlqq + m\L*L + mYti+ ^{H* Hf - As {HQ*q + H*Qq*) ~ Al [HL*1 + H*LI* 



9i 



9l 



{h*HQ*Q + H*Hqq + Q*Qqqj + {h*HL*L + H* Hl*l + L*Ll*lj + ^ ( H*H - Q*Q 
(h*H- + ^ (^H*H + Q*Q -2q*qy + ^ (^H* H + L* L - 2l*t ^ 



(52) 



where Q and q are squarks, L and £ are sleptons, and 
iJ is a Higgs boson, although again for simplicity, we 
take these fields to be real. One local minimum of this 
potential is at (^Q^ = (q) = (^Lj = (^tj = and (H) = 

nih/V^- For these values of the coupling constants: 



Potential Along Line Connecting Minima 



WV2 GeV 



itLq = niq = 15 V 2 GcV 
As = 31 GeV 
91 = 92 = -6 



niH = .5\/2 GeV 
Al = 23 GeV 
A = .006 

y = 1, (53) 



the global minimum is at the minimum mentioned above, 
while a local minimum occurs at (Q) = (q) = 0, (L) = 



5.025 GeV, {Ij = 5.136 GeV, and (H) = 8.292 GeV. 

Since the squarks do not acquire a VEV in this vacuum, 
baryon number is conserved and can be used to construct 
Q-balls. 

In the false vacuum, the lightest squark q' has a mass 
rriq' of 7.80 GeV. Because lepton number is not conserved, 
the sleptons mix with the Higgs boson; the lightest of 
these eigenstates is h' with a mass of 6.83 GeV. We as- 
sume that the quark acquires a mass from the term yqHq 
in the Lagrangian; then it has mass y (H) = 8.29 GeV, 
and so the lightest squark is stable against decay into a 
quark. 

The potential along the line connecting the minima is 

f7((/)) = .32^^ - 7.38 GeV 0^ + 41.2 GeV^ c/)^, (54) 

which is shown in Fig. [Tj We notice that the barrier 
separating the local minimum from the global minimum 
is quite large, which dramatically suppresses tunneling 
through thermal fluctuations. By expanding the poten- 
tial in terms of the appropriate eigenstates in the false 
vacuum, one can show that Q-balls constructed of q' and 



> 



300 



200 - 




3 14 



FIG. 7. The potential along the line connecting the false 
vacuum to the true vacuum. 



h' fields do exist. Thus we have all of the necessary in- 
gredients for a solitosynthesis-induced phase transtion. 

While most of our analysis could be straightforwardly 
applied to this situation, the Bethe-Salpeter regime anal- 
ysis cannot. Firstly, because of the balance between the 
squark and quark masses, we do not have ruh' ^ rriq' 
as required for the Bethe-Salpeter equation. As yet, the 
Bethe-Salpeter equation for the exchange of massive par- 
ticles is unsolved. 

Secondly, as was noted in the subsection "Bethe- 
Salpeter Regime" of the section "Properties of Q-balls 
in the False Vacuum", the field content of states de- 
scribed by the Bethe-Salpeter equation does not neces- 
sarily match the field content of Q-balls described by the 
thick wall approximation. In our numerical example, we 
made the difference between the two small by choosing 
the global minimum such that sin(0) w 1. However, in 
this scenario, we must have sin(f?) <^ 1 because we must 
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tunnel to a state near the global minimum, which has 
{q') = 0. Therefore, the analysis should be modified to 
account for the fact that the Q-balls consist almost en- 
tirely of the h' field, with very little of the q' field. 

However, even if such a phase transition could not have 
occured in the evolution of our universe, it is still impor- 
tant to study the regions of parameter space in the MSSM 
in which such a phase transition could occur. Such a 
phase transition could destabilize a vacuum previously 
thought to be stable on cosmological timescales, leading 
to further constraints beyond those of [23] . 



CONCLUSIONS 

We have considered the properties and growth of Q- 
balls in the false vacuum in each of four regimes, rang- 
ing from extremely small Q-balls to extremely large Q- 
balls, and we have demonstrated that phase transitions 
induced by solitosynthesis are indeed possible. While we 
have used a toy model inspired by the MSSM, such phase 
transitions occur what any model with a similar potential 
and vacuum structure. 

These phase transitions are of cosmological interest. 
One unique aspect of these phase transitions is that the 
resulting vacuum carries a net charge afterwards; as men- 
tioned in the introduction, this has been suggested as 
a possible baryogenesis mechanism in the AfHeck-Dine 
mechanism. However, more work must be done to es- 
tablish that these phase transitions occur in the most 
promising MSSM potentials for this scenario. Even if 
this cannot be established, the existence of such phase 
transitions may set new bounds on the allowed MSSM 
parameter space. 

The author would like to thank Alex Kusenko for very 
helpful discussions. This work was supported in part by 
DOE grant DE-FG03-91ER40662. 



Appendix: The Bethe-Salpeter Equation 

In this appendix, we will review the Bethe-Salpeter 
equation and derive the results given in the text for un- 
equal masses in more detail. The Bethe-Salpeter equa- 
tion describes the relativistic bound states of a strongly 
interacting system, described by the ket \B). The am- 
plitude is: 



(A.55) 



where P is the four-monentum of the bound state. One 
may simplify this equation by the ladder approximation, 
in which one considers only diagrams of the type shown in 
Fig. [8] The Bethe-Salpeter equation has been solved in 
the Wick-Cutkosky model, in which the particles involved 
are two scalars exchanging massless quanta [IS;. In this 



I I I 



I I I 



FIG. 8. These Feynman diagrams are considered in the ladder 
approximation of the Bethe-Salpeter equation. 



model, the Bethe-Salpeter equation after a Wick rotation 
can be written as: 



M 



M 

Y " 

167r4 



[p - ky + ' 

(A.56) 



where M is the mass of the bound state, m is the mass 
of the (j) fields, and A is the coupling between (j) and the 
massless particle. The momenta of the particles on the 
top and bottom are P/2 ± p. The wavefunction can be 
found by solving a Friedholm differential equation pi] ; 
however, we need only the bound state masses, which are 
given by [H]: 



a 



2 \ 



M„ = 2m(l-^J, 



(A.57) 



if a A^/lQirm^ < 1. 

In the Wick-Cutkosky model, the particles at the top 
and bottom of the ladder are identical; however, this can 
be weakened. Maintaining the ladder approximation, the 
Bethe-Salpeter equation for non-identical scalars inter- 
acting through massless exchange is: 



[(m + A)2 + {p- i77iP)2] [{m - Af + {p + i7]2Pf] Mp) 

^' ^"^'v^, (A.58) 



(p - kf 



where the masses of the particles on the top and bottom 
are to ± A. The coupling constant A' is related to the 
coupling constants at the top and bottom interactions by 
-A' = gtopgboUom/^(^T^^ ■ Vi ^-nd ?72 come from transform- 
ing to a "center of momentum" reference frame; these 
are: 



nit 



op 



om 



^bottom 



nitop + nibottom 



The energy-momentum four-vector of the bound state, 
P, is given by (0, M) where M is the bound state mass. 

As derived in [12], this equation can be related to the 
Wick-Cutkosky model Bethe-Salpeter equation. The re- 
sult is that if A' = F{M'^) in the case of equal masses. 
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then for unequal masses: 
A' 

F 



1 - Is? I'm? 



M2 - 4A2 

1 - A2/to2 



(A.59) 



We will use this to derive the ground state for unequal 
masses. First, we note that for equal masses A' = A^ / 
since gtop — gbottom — A. Therefore, equation 
(|X57| is: 



2to ^ ^ 8m4 



(A.60) 



Squaring this and keeping the lowest order terms using 
Al <^ m? gives: 



Am? 

which can be solved for: 



= 1 - 



4m4 



A' = 



1 - 



4m" 



(A.61) 



(A.62) 



Then for unequal masses: 



(1- A2/m2 



'4m4 1 AP - 4A2 



A 



12 



4m4 



and solving for M'^ gives 



1 - 



4m2 1 _ /S?/n? 



1 M2 - 4A2 



M2 ^ 4A2+4m2 1 



A2 

rn? 



Am? 1 - A2/r 



1 - 



(A.63) 
(A.64) 



1 



4m4 (l-A2/m2)2 
(A.65) 



as in (26). Setting A = restores the equal mass case. 
This remains valid as long as A' <^ m^ . If this fails, 
however, we are in a strong coupling regime in which 
the Bethe-Salpeter equation underestimates the binding 
energy; therefore Q-balls will be more likely to bind to- 
gether than we find here. Thus if we find a phase transi- 
tion with this approximation, the phenomenon will still 
occur if the binding energies were calculated exactly. 
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